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ABSTRACT 


Efficient underwater acoustic communication and target locating systems 
require detailed study of acoustic wave propagation in the sea. Many 
investigators have studied the absorption of acoustic waves in ocean water 
and formulated empirical equations such as Thorp’s formula, Schulkin and 


Marsh model and Fisher and Simmons formula. The Fisher and Simmons 
formula found the effect associated with the relaxation of boric acid on 
absorption and provided a more detailed form of absorption coefficient 
which varies with frequency. However, no simulation model has made for 
the underwater acoustic propagation using these models. This paper reports 
the comparative study of acoustic wave absorption carried out by means of 
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Absorption modeling in MATLAB. The results of simulation have been evaluated using 
Boric acid measured data collected at Desaru beach on the eastern shore of Johor in 
Frequency Malaysia. The model has been used to determine sound absorption for given 
values of depth (D), salinity (S), temperature (T), pH, and acoustic wave 
transmitter frequency (f). From the results a suitable range, depth and 
frequency can be found to obtain best propagation link with low absorption 

loss. 
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1. INTRODUCTION 

Increased interest in defense applications, off-shore oil industry, and other commercial operations 
provides a motivation for research in signal processing for the underwater environment. In the underwater 
environment, acoustics waves are more practical for applications such as navigation, communication, and 
other wireless applications due to the high attenuation rate of electromagnetic waves. Acoustic propagation is 
characterized by three major factors: attenuation that increases with signal frequency, time-varying multipath 
propagation, and low speed of sound (1500 m/s) [1]. 

The fading in underwater environment depends on the distance, frequency and sensors location, and 
can be divided into long term fading and short term fading. Long term fading caused by sound propagation is 
effected by spreading loss and absorption loss which is a function of range and frequency. Short term fading 
that composed of multipath and Doppler spread which is a random function of distance and time. 

As the attenuation of sound in the ocean is frequency dependent, the ocean acts as a low-pass filter 
for ambient noise and the underwater systems operate at low frequencies, for example, on the order of tens of 
kHz [1]. No two deployment regions within the ocean with have the same depths ranging from tens of meters 
to a few kilometers with node placement that varies from one network to another [2]. 
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Acoustic energy can be reflected and scattered from both the surface and bottom of the ocean as 
shown in Figure 1., permitting the possibility of more than one transmission path for a single signal [1, 2]. 
The phenomenon is known as multipath fading. The overlapping of consecutive pulses in digital 
communication results in intersymbol interference (ISI) that increases the bit-error rate in the received data 
[3-5]. Therefore, underwater data communication links generally support low data rates [1, 3]. This study 
aims to show the effects of absorption loss in propagation link between the transmitter and receiver. 
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Figure 1. Underwater acoustic environment [1] 


2. RESEARCH METHOD 
2.1. Sound Speed 

The speed of sound in seawater is a fundamental oceanographic variable that determines the 
behavior of sound propagation in the ocean. Many empirical formulas have been developed over the years for 
calculating sound speed using values of water temperature, salinity, and pressure/depth. A simplified 
expression for the sound speed was given by Medwin [6]: 


c = 1449.2 + 4.6T — 0.0557? + 0.00029T? + (1.34 — 0.01T)(S — 35) + 0.016d (1) 


where c is the speed of sound in seawater, T is the water temperature (in degrees Celsius), S is the salinity (in 
parts per thousand) and d is the depth (in meters). 


2.2. Total path loss 

An acoustic signal underwater experiences attenuation due to spreading and absorption. Path loss is 
the measure of the lost signal intensity from the projector to the hydrophone. Spreading loss is due to the 
expanding area that the sound signal encompasses as it geometrically spreads outward from the source [7]. 


Pl spreading (R) =k * 10log(R) dB (2) 


where R is the range in meters and k is the spreading factor. When the medium in which signal transmission 
occurs is unbounded, the spreading is spherical and the spreading factor k = 2; whereas in bounded 
spreading, it is considered as cylindrical k = 1. In practice, a spreading factor of k = 1.5 is often considered 
[2]. 

The absorption loss is a representation of the energy loss in form of heat due to viscous friction and 
ionic relaxation that occur as the wave generated by an acoustic signal propagates outwards; this loss varies 
linearly with range as follows [7]: 


Prod f) = 10 log(a(f)) x R dB (3) 


where r is range in kilometres and a@ (f) is the absorption coefficient. The absorption coefficient using 
different emperical models will be discussed in next section. 
Total path loss is the combined contribution of both the spreading and absorption losses [1]. 
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Total Path Loss = k * 10 log(R) + 10 log(a(f)) * R (dB) (4) 


2.2. Absorption Loss Empirical Models 
The acoustic energy of a sound wave propagating in the ocean is partly: 


a. absorbed, 1.e. the energy is transformed into heat. 
b. lost due to sound scattering by inhomogeneities. 


On the basis of extensive laboratory and field experiments the following empirical formulae for 
attenuation coefficient in sea water have been derived. There are three formulaes: Thorp’s formula, Schulkin- 
Marsh model and Fisher-Simmons formula and the frequency band for each one as shown in Figure 2. 

R. E. Francois and G. R. Garrison [11] have formulated the equation for the sound absorption in the 
frequency range 400Hz to 1MHz which includes the contribution of Boric Acid, Magnesium Sulfate and Pure 
water. The results given by this equation are very close to the practical results. 
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Figure 2. Diagram indicating empirical formulae for different frequency domains 


The absorption coefficient for frequency range 100 Hz to 3 kHz can be expressed empirically using 
Thorp’s formula [2, 8] which defines a [dB/m] as a function of f [kHz] 


= P A —4¢2 =3 
a(f) = (0.11 + 4442.75. 10-4? + 0.0003) 10 (5) 


The absorption coefficient for frequency range 3 kHz to 500 kHz can be expressed empirically using 
Schulkin and Marsh model [2,8] which defines a [dB/m] as a function of f [kHz] 


Afr f4 Bf 
2 — ie 
ftf fp 





a(f) = 8.686. 103 ) (1 — 6.54. 1074P). 1073 (6) 


where A and B are constants A = 2.34.107”*,B = 3.38.107°, P is the hydrostatic pressure [kg/cm3], 


1520 


fr = 21.9. 10!6- (F273)! [kHz] is the relaxation frequency. 
An alternative expression for the absorption coefficient a(f) [dB/m] is given by the Fisher and 
Simmons formula in range of frequency 10 kHz- 1 MHz [9] 


(2 (2 
aA =| AP 7h + At zp 


o] D 
B(OH)3,Boric acid MgSOy¿,magnesium sulphate 


f + AR f?  |.1073 (7) 


H2O,pure water 


where A}, Az, and Azare the coefficients represent the effects of temperature, while theP,, P,and 
Pz coefficients represent ocean depth (pressure) and f,, fọ represent the relaxation frequencies of Boric acid 
and (MgSO4) molecules [9]. As a result, the absorption coefficient is proportional to the operating frequency. 
1. Boric acid B(OH); 


8.686 
A, =—— 


109:78Ph-s 
P, = 1 
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2.8 (101273 
= Z. — T +273 
fi 35 


2. Magnesium Sulphate MgSO4 
S 
A, = 21.44.—. (1 + 0.025T) 


Py = 1 — 1.37. 10~4 Zmar + 6.2. 107°. ZZ ax 
g—__1990 
_ 8.17. (10 @-273)) 
* 1+ 0.0018(S — 35) 


3. Pure water H,0 


_ ( 4.937.10-4* — 2.59.10-5T + 9.11.1077T? — 1.5.10-8T? for T < 20°C 
3 13.964. 1074 — 1.146.107 + 1.45.1077T2 — 6.5.1078T3 for T > 20°C 


P3 = 1—3.83.10~5Zmax + 4.9. 107°. Z2 ax 


with f in [kHz], T in [°C], S in [ppt]. And where Z,,,,, pH and c denote the depth in [m], the pH-value and 
the sound speed in [m/s] respectively. 

R. E. Francois and G. R. Garrison [11] have formulated the equation for the sound absorption in the 
frequency range 400Hz to 1MHz which includes the contribution of Boric Acid, Magnesium Sulfate and Pure 
water. The results given by this equation are very close to the practical results. 
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Figure 3. Experiment test site 


A general diagram showing the variation of a(f) using Francois and G. R. Garrison formula with the 
three regions of Boric acid, B(OH)3, Magnesium sulphate, MgSO, and Pure water, H2O is depicted in Figure 
4. It can be observed that for the Boric acid region, Attenuation is proportional to f * . And for the regions 
Magnesium sulphate and pure water also Attenuation is proportional to f ° . In the transition domains it is 
proportional to f. 

With fixed the range at 1000 m, the path loss is proportional to the operating frequency and it is 
increased when the operating frequency increased and the contribution of the absorption term is less 
significant than the spreading term as shown in Figure 5. 

Figure 6. and Figure 7. with fixed the range at 10Km and 100 Km respectively. As range increases 
and the absorption term begins to dominate, any variations in œ also becomes more significant. For data 
communication, the changes in the attenuation due to signal frequency are particularly important as the use of 
higher frequencies will potentially provide higher data rates. 

The total path loss versus the range at different frequencies shown in Figure 8. The path losses 
Increase significantly at high frequency compared with low frequency. Therefore, the attenuation of sound in 
the sea is frequency dependent, the sea acts as a low-pass filter. 
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Figure 6. Total path losses as a function of frequency with range 10 Km 
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Frequency VS. Total path losses 
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Figure7. Total path losses as a function of frequency with range 10 Km 
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Figure 8. Total path losses as a function of range 


3. CONCLUSION 

In this paper, the losses due to phenomena of acoustic wave propagation and absorption have been 
studied. The absorption coefficients in dB/km vs signal frequency for Francois and G. R. Garrison formula 
shows that in general @ increases with increasing frequency at any fixed temperature and depth. For very 
short-range communication, the contribution of the absorption term is less significant than the spreading 
term. As range increases, the absorption term begins to dominate, any variations in æ also becomes more 
significant. For data communication, the changes in the attenuation due to signal frequency are particularly 
important as the use of higher frequencies will potentially provide higher data rates. Nevertheless, no other 
type of radiation can compete with low-frequency sound waves to communicate over great distances in 
oceanic environment. 
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